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We show, by using Regge calculus, that the entropy of any finite part of a Rindler horizon is, in 
the semi-classical limit, one quarter of the area of that part. We argue that this result implies that 
the entropy associated with any horizon of spacetime is, in the semi-classical limit, one quarter of 
its area. As an example, we derive the Bekenstein-Hawking entropy law for the Schwarzschild black 
hole. 
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' A remarkable discovery was made by Bekenstein and Hawking almost 30 years ago when they found that black 
' hole has a certain entropy which is one quarter of its horizon area[Q ||. As a consequence, black holes emit thermal 
radiation with a spectrum which is essentially that of a black body. After Bekenstein and Hawking it was soon 
I discovered by Unruh||, ^, || that not only black holes but also the so-called Rindler horizon of an accelerated observer 
emits thermal radiation. So it appeared that the black hole is not necessary for the creation - from the point of view 
^ ' of a certain observer - of thermal radiation out of the vacuum. Rather, what is essential for the creation of thermal 
, radiation, is the existence of a horizon. In short, it appears that all horizons of spacetime, no matter whether they 
are black hole, Rindler, de Sitter, or cosmological horizons, have certain thermodynamical properties in common. For 
, instance, all horizons of spacetime emit thermal radiation, fn particular, it appears - and this is one of the main 
■ theses of this paper - that if we consider any part of any spacetime horizon, then one can associate with that part an 
I I entropy which is one quarter of the area of the part. 
^ . In this paper we give support to this thesis by means of a detailed analysis of the thermodynamical properties 
00 of Rindler horizons. A novel feature in our analysis is the use of the Regge calculus^ |^, ^ approach to Einstein's 
CN general relativity. In short, Regge calculus is an approach to general relativity, where spacetime is modelled by a 
• piecewise flat, or simplicial, manifold and the geometrial properties of that manifold are coded into its edge lengths. 
1^ ' Using Regge calculus we show that in the semi-classical limit each part of the Rindler horizon of spacetime possesses 
an entropy which is one quarter of its area. Finally, at the end of this paper, we shall show how the analysis used in 
the study of Rindler horizons may also be used for the derivation of the Bekenstein-Hawking law for the black hole 
entropy. 

I Rindler horizons are such horizons of spacetime which appear in the rest frame of an accelerated observer. In 
I , general, the world line of a uniformly accelerated observer in flat two-dimensional Minkowski spacetime has the 
^Jpj' equation (unless otherwise stated we shall always have c = G = fi — ~ 1)Q 



o 
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X 



X2-T2 = 1, (1) 



a 



JH \ where a is the proper acceleration of the observer, and X and T, respectively, are the niinkowskian space and time 
coordinates. The world line of the observer may be written in the parametrized form 

T(i) = isinh(ai), (2) 
a 



= -cosh(at). (3) 
a 

In this expression, t is the proper time of the observer. The world lines of the two observers, being uniformly accelerated 
in opposite directions, have been drawn in Fig. |l]. In that fi gure we may also see the Rindler horizon (which, in the 
observer's rest frame is, in our two-dimensional figure, represented by a point) of the accelerated observers. 

Consider now two oppositely accelerated observers. The spacelike hypersurfaces where the proper times t of the 
observers are ti and t2 (ti < are orthogonal to the world lines of the accelerated observers. So far we have just 
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FIG. 1: Rindler spacetime. The hypersurfaces where t — constant are orthogonal to the world lines of the accelerated observers. 
The curvature of these hypersurfaces is concentrated at the Rindler horizon at time T — 0. 



recalled the standard properties of accelerated observers in flat spacetime. At this point we turn to quantization. We 
evaluate the propagator 

K{qab{t2),t2;qab{tl),ti). 

This propagator gives the propability amplitude that the metric tensor on the hypersurface t — t2 is qah{t2) provided 
that the metric tensor on the hypersurface t — ti was qab{ti)- Formally, the propagator may be written as a path 
integral: 

K{qab{t2)it2\qab{ti),h) = J '^ig^.v] exp {iScig^^]) , (4) 

where the integration has been performed over the four-metrics g^i, between the hypersurfaces t ~ ti and t ~ t2 
(shaded region in Fig. ^) . Sq [gfj.i^] , in turn, is the classical gravitational action: 

SG[gtiA = T^l RV^d^x+SdM, (5) 

where R is the Riemannian scalar and g is the determinant of the four- metric g^,^. In the first term on the right hand 
side of Eq. the integration is performed over the spacetime region between the surfaces t — ti and t — t2, whereas 
the second term is the boundary term to the gravitational action. 

At the present state of research it is not known how to evaluate explicitly the right hand side of Eq. (^. However, it is 
possible to approximate the propagator by means of the so-called semi-classical approximation. In this approximation 
we may write the propagator as 

K{qab{t2),t2;qab{tl),ti) W [(Jafc 1 ) , ^afc (^2 )] exp (iS'g [.g^tJ.] ) , (6) 

where N^qab{ti), qab{t2)\ is some slowly-varying functional of qab{ti) Siud 906(^2)1 a-nd Sq is the gravitational action 
when Einstein's field equations are satisfied. For flat spacetime R = 0, and therefore we are left with the boundary 
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FIG. 2: The hypersurfaces t = ti and t = t2 bound a certain region of Rindler spacetime (shaded region in this figure). When 
calculating the propagator we integrate over all four-metrics g^i^ in this region. 



terms only. It is possible to choose the boundary term at asymptotic infinity such that it vanishes in flat spacetime, 
and therefore all contribution to the boundary terms comes from the hypersurfaces t = ti and t = t2. It is at this 
point where Regge calculus enters the stage. 

In Regge calculus the gravitational action is written as 

3 k 

In this expression, Aj's are the areas of those two-simplices, or triangles j, which do not lie on the boundary of 
spacetime, and e^'s are the corresponding deficit angles. A^'s, in turn, are the areas of those triangles k which 
lie on the boundary of spacetime, 0fc's being the corresponding deficit angles. If the triangles k lie on a spacelike 
hypersurface of spacetime, then 0^ may be understood as the boost angle between the timelike normals of the two 
tetrahedra having the triangle k in common. If spacetime is flat, we have 

= (8) 

for every j, and the gravitational action reduces to 

k 

In other words, if we have a spacelike hypersurface embedded in flat spacetime such that its curvature is concentrated 
along certain flat two-spaces, we get the action - up to the term (Stt)"^ - by just multiplying the areas by the 
corresponding boost angles, and summing the products. This result is particularly useful in our present study of 
Rindler horizons: Rindler horizons are flat two-spaces, and therefore they may be constructed from flat triangles. 
Morever, as one can see from Fig. ^ the curvature on the hypersurfaces t = constant is concentrated along the 
Rindler horizon, being zero elsewhere. Even more interesting is the fact that the boost angles are the same for all 
triangles on the Rindler horizon. Therefore the classical action Sg takes, when spacetime is flat, the form 

Sg = ^A(</>2 - 0l), (10) 
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where A is the total area of the considered part of the Rindler horizon, whereas and 412, respectively, are the boost 
angles between the future pointing tangent vectors of the world lines of the oppositely accelerated observers at times 
ti and t2. Therefore we may approximate the propagator as 

K{qab(t2)>h;qab{ti),ti) » N[qab{tl),qab(t2)]'^^P 

We are now prepared to enter the thermodynamics of Rindler horizons. As it is well known, thermodynamics is 
just quantum mechanics in euclidean spacetimej^. More precisely, if a system is periodic in euclidean spacetime with 
period P such that after the elapsed euclidean time /3 it returns to the original point x in configuration space, then 
the partition function of the system is 

Z{P) = J dxK{x,-iP;x,0). (12) 

In other words, we have replaced the lorentzian time t in the propagator by the term —if3, where /3 is the euclidean 
time. Thermodynamically, the quantity (3 may be intepreted as the inverse temperature of the system. We shall use 
this idea when we calculate the partition function, and thereby the entropy, of the Rindler horizon. 

To begin with, consider an accelerated observer in euclidean spacetime. We just replace in Eqs. and (j3|) T by 
—iT' and t by — ir, where T' and r are euclidean time coordinates. We get: 

T' = isin(aT), (13) 
a 



Hit 



A{(t>2 - (1)1 



(11) 



X ^ — cos(aT). 
a 



(14) 



In other words, the world line of an accelerated observer in lorentzian spacetime is, in euclidean spacetime, a circle. 
In Fig. H we have drawn the world lines of two accelerated observers in euclidean spacetime. The hypersurfaces t = 
constant are replaced by the hypersurfaces where r = constant. It is remarkable that spacetime geometry is now 
periodic with respect to the euclidean time r: An observer returns to the same point in euclidean spacetime after the 
elapsed euclidean time 



a 



(15) 



and therefore the temperature experienced by the accelerated observer is exactly the Unruh temperature of the Rindler 
horizon iSl: 



~ ^ - 2^- 



(16) 



In contrast to spacetime with lorentzian signature, the hypersurfaces r — constant sweep over the whole spacetime. 
However, during just one period they sweep the spacetime two times. Because of that we must divide the action by 
two when we calculate the partition function of the Rindler horizon from the point of view of a single observer. 

We now proceed to calculate the partition function. To begin with, we consider the action in euclidean spacetime. 
In euclidean spacetime we must replace the boost angles 0i and (j)2, respectively, by the quantities —iOi and —i02, 
where 61 and 62 are just the "ordinary" euclidean angles between the normals of the tetrahedra meeting at their 
common two-simplex on the Rindler horizon. Hence, the euclidean propagator takes, in the semi-classical limit, the 
form: 



KE(qab{r2),T2;qab{Ti),Ti) fa N[qab{Ti),qab{T2)]exp - 77-^(^1 ~ (^2) 



(17) 



During one period the angle between the normals of the tetrahedra meeting at this common two-simplex decreases 
by 47r, and therefore, taking into account that the action must be divided by two for the reasons explained before, 
the partition function is, in the semi-classical limit. 



V 



^ab^^) N qab(^^yqab{0) CXp jAj. 



'2n 



1 



(18) 



FIG. 3: In euclidean spacetime the world line of a uniformly accelerated observer is a circle. The hypersurfaces where the 
euclidean proper time r of the observer is constant sweep over the whole spacetime. 



Since gab(O) 



Qabi^), we integrate over a single variable qab{^), and when this variable is integrated out, we get 




where A/" is a constant. In other words, we have obtained an exciting result that, at least in this semi-classical limit, 
the partition function of the considered part of the Rindler horizon has the negative of one quarter of the area of that 
part in the exponential. 

To calculate the entropy of the Rindler horizon we must be able to express the partition function in terms of the 
euclidean period /?. To this end, consider what an accelerated observer actually observes about the properties of 
his Rindler horizon. Such an observer sees that, in his frame of reference, all bodies tend to fall, with a uniform 
acceleration o, towards the Rindler horizon. According to the Principle of Equivalence the observer has no means 
to decide whether he is in a uniformly accelerating motion, or in a uniform gravitational field caused by a mass 
distributed uniformly along a plane. According to Gauss' law for gravitational fields the mass enclosed by the closed 
two-surface a is given by a surface integral 



where g is the gravitational field. If the observer is accelerated along the a;-axis, he may conclude that he is in the 
gravitational field 



and therefore, because the Rindler horizon lies in the yz-pl&ne, the mass contained by the part of area A of the Rindler 
horizon, is 




(20) 



ai, 



(21) 



M 




(22) 
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which follows from a consideration similar to the one performed in the standard textbooks of electromagnetism, when 
the charge density of a plane is calculated by means of Gauss' law for electric fields. In other words, the area of the 
considered part of the horizon may be written in terms of its observed mass M and inverse temperature /3: 

A = 2(3M, (23) 

and therefore the partition function Z may be written, in the semi-classical limit, in terms of the average value (M) 
of the mass: 

Z = AAcxp(-i/3(M>). (24) 
However, the average value of M may calculated from the partition function: 

{M)=-^\nZ. (25) 



This gives us a differential equation for (M): 



which, in turn, has the general solution 



= (M). (26) 



(1/) = C/3, (27) 
where C is a constant. Hence, the partition function takes, in terms of /3, the form 

Z(/3)=AAexp(-ic/32). (28) 

We now proceed to calculate the entropy of the Rindler horizon. In general, the entropy of any system may be 
calculated from the partition function: 

S = -p(^^lnz)+lnZ. (29) 



This gives us (we take Af ~ 1): 



and it follows from Eqs. (E3|) and (|27|) that 



S = ^C(3\ (30) 



S=Ia. (31) 



In other words, we have obtained a result which states that, from the point of view of an accelerated observer, the 
entropy of any finite part of the Rindler horizon is one quarter of the area of that part. As such our result is similar 
to the Bekenstein-Hawking law for black hole entropy - the only difference is that the area of the event horizon of the 
black hole has been replaced by the area of the considered part of the Rindler horizon. 

Actually, our result may be used to show that any horizon of spacetime has an entropy which is one quarter of its 
area. That is because, locally, any small enough part of any horizon may be considered, in an appropriate coordinate 
system, as a part of a Rindler horizon. The entropy of the horizon as a whole, in turn, is the sum of the entropies of 
its parts, and since the horizon may be constructed from tiny Rindler horizons, each having an entropy equal to one 
quarter of its area, the entropy of the whole horizon is one quarter of its total area. 

To see what all this means consider, as an example, the Schwarzschild horizon. The Schwarzschild metric in 
Schwarzschild coordinates t and r is 



ds^ ^-{l- '^)df + + r\de^ + sin2 0d</)2), (32) 
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where M is the mass of the Schwarzschild black hole. An observer at rest with respect to the coordinates r, 9 and 4> 
observes that the hole radiates thermal radiation with a characteristic temperature Q 

Th =: (33) 



where we have taken into account the "red shift" of the temperature. In other words, an observer at rest very close 
to the horizon r = 2M may measure an infinite temperature for the black hole radiation. 

Consider now the radiation process from the point of view of an observer in free fall, very close to the horizon. A 
propriate coordinates for the study of the properties of freely falling observers are the so-called Novikov coordinates 
(for a detailed survey of the Novikov coordinates, see Appendix). In terms of the Novikov coordinates (r, i?*) the 
Schwarzschild line element takes the form 

ds' = -dr' + + ,2 (^^2 ^ ^^^2)_ (34) 

For an observer in radial free fall such that ^ = when t = 0, R* is constant and r is the proper time of such an 
observer. At the bifurcation point in the Kruskal diagram of Kruskal spacetime R* — t — 0, and for an observer in 
radial free fall through the bifurcation two-surface R* ~ 0. Close to the bifurcation two-surface we may write as 
a Taylor expansion (see Appendix): 

^ = 4Mi?* + 3Mi?*^ + + . . . , (35) 

and therefore the spacetime metric becomes, at the bifurcation two-surface, 

ds = -dr' + [AMfdR*' + {2Mfdn' (36) 

where d^' is the line element on a unit two-sphere. If we rescale the Novikov coordinate R* such that we define a 
new coordinate 

R*:=AMR\ (37) 

we get 

ds' = -dr' + dR*' + {2Mfdn'. (38) 

In other words, the Novikov coordinate system provides, at the bifurcation point, a coordinate system which is locally 
flat, or locally inertial. 

It is instructive to compare the expression ( ^8| ) to the metric given by the Kruskal coordinates u and v. In general, 
the Schwarzschild line element takes, in Kruskal coordinates, the form 



ds' _ ^^ 

and at the bifurcation two-surface we have 



/32M^ 



y--/^^'\-dv' + du') + r'da', (39) 



ds' = ^^{-dv' + du') + {2M)'dn'. (40) 
e 

Hence we observe that very close to the bifurcation two-surface between the Kruskal coordinates (u, u) and the 
coordinates {t,R*) there is the relationship 

r^i^.. (41) 



R* = ^u. (42) 

In other words, close to the bifurcation two-surface the coordinates r and R* are, up to the rescaling factor AM/y/e, 
identical to the Kruskal coordinates v and u. 
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We are now coming to the similarities between the Rindler and Schwarzschild horizons. As it is well known, the 
relationship between the Schwarzschild and Kruskal coordinates is, when r > 2M0: 

,= (^^lY\rm^^^(t\ (43) 

V2M J \AMJ' ^ ' 

V2A/ / V4M/ ^ ' 
To consider this transformation very close to the bifurcation two-surface, we define 



r 

and the transformations of Eqs. (Esl) and (H) become, for small 8 



T = AM (5 sinh 

^*=4M5cosh(^j, (47) 
where we have defined a new time coordinate 

d := tb, (48) 

which actually is the pro per time of an observer at rest with respect to the Schwarzschild coordinates. Comparing 
Eqs. (^) and (^) to Eqs. (^) and (^) we notice that observers at rest with respect to the Schwarzschild coordinates 
are, quite simply, observers which are accelerated with respect to the observers in free fall with a proper acceleration 

(49) 



4M<5 4M,/rnM 



In other words, the appearance of the Schwarzschild horizon from the point of view of an observer at rest with respect 
to the hole is simply due to the fact that such an observer is accelerated with respect to a locally inertial observer. In 
this sense there is no difference between a Rindler horizon and any sufficiently small part of a Schwarzschild horizon. 
Both appear simply because the rest frame of the observer is non-inertial. From the point of view of an inertial 
observer very close to the horizon, no horizon, nor any of its effects, such as Hawking radiation, do appear. 

From this point on, the thermodynamical properties of the Schwarzschild horizon may be obtained in exactly the 
same way as we obtained the thermodynamical properties of the Rindler horizon. For instance, the temperature of 
the horizon is 

Th = ^ = (50) 



which is Eq. (pS]). Moreover, any sufficiently small part of the horizon has an entropy which is one quarter of its area, 
and therefore, a jortirior% the entropy of the hole is one quarter of its horizon area. A similar chain of reasoning may 
be applied to any horizon of spacetime. 

In this paper we have considered the thermodynamical properties of spacetime horizons in general. Our consid- 
erations were based on detailed investigations of Rindler horizons. We argued, by using Regge calculus that, in the 
semi-classical limit, one may associate with any finite part of a Rindler horizon an entropy which is one quarter of 
the area of that part. This is the entropy observed by an accelerated observer. Moreover, we showed that the results 
originally obtained for Rindler horizons may be applied for the derivation of the thermodynamical properties of the 
Schwarzschild horizon as well. This conclusion was based on the observation that an observer at rest with respect to 
the horizon is accelerated with respect to an observer in free fall through the horizon. In other words, any sufficiently 
small part of a Schwarzschild horizon may be considered as a part of a Rindler horizon. Since the entropy of the 
Rindler horizon is one quarter of its area, so is the entropy of the Schwarzschild horizon as well. A similar chain of 
reasoning may be applied to any horizon, and therefore it appears that one may argue, in the semi-classical limit, that 
the entropy of any spacetime horizon is, from the point of view of an observer at rest with respect to that horizon, 
one quarter of its area. 

Note: During the preparation of this paper we became aware of an interesting recent paper written by 
Padmanabhanlin]. He also showed that the entropy of any horizon is one quarter of its area. His methods, however, 
are completely different from ours. 
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APPENDIX A: NOVIKOV COORDINATES CLOSE TO THE BIFURCATION POINT 

1. Novikov Coordinates 



When one is considering observers in free fall towards the Schwarzschild black hole, it is often useful to apply 
Novikov coordinates. In this coordinate system a freely- falling observer at rest with respect to the Schwarzschild 
coordinates, for the Schwarzschild time t = 0, remains at rest during his whole journey. At the time t = the 
observer is released in radial free fall and the Schwarzschild radial coordinate r of the observer can be written at the 
time t = as[0 

r = r^a.^2M{l + R*^), (Al) 

where R* is the radial coordinate in Novikov coordinates. The time coordinate in Novikov coordinates is the proper 
time T of this freely-falling observer. Thus, in the Novikov coordinate system spacetime points are identified by means 
of the radial coordinate R* and the proper time coordinate r. 

The equation of motion for an observer in free fall in Schwarzschild spacetime is 



1 = 



C2 



,-,2 



1_ 2M ^_ 

r I 

where C is a constant and the dot stands for the proper time derivative. Hence we get 



2M 



r 



It is easy to see that the observer is at rest with respect to the Schwarzschild coordinate r when 

2M 



C2 - 1' 



From equations (Al) and (A3) we may solve f in terms of R* and r: 



± 



2M 

r 



1 



-I- 1 

Furthermore, we can express the proper time t as an integral 



1/2 



± 



2M 



1 



R* 



1 



-1/2 



dr , 



which gives the relationship between the Schwarzschild coordinate r and the Novikov coordinates R* and r: 

(l + i?*2) 



± 



2M 



2M AK'P{R*'^ + l 
+ {R*'^ + 1)3/2 j^rccos 



1/2-1 



2M{R 



*2 , 



1 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



The Schwarzschild line element can also be written in terms of the Novikov coordinates R* and r, and it takes the 
form 

+ 1 / dr^^ 



-dr' 



R 



*2 



dR* 



dR*^ +r\de^ +s\n^6d(l)^). 



(A8) 
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Here r should be thought as a function r(T, R*) given imphcitly by Eq. ( A7). The point r = i?* = is the bifurcation 
point u = V = of Kruskal spaceti me. A t this point r = 2M. It is possible to obtain the Taylor expansion for (^^fir) 
and thereby write the line element (A8) close to the bifurcation point in terms of R* and r. 



Let us write Eq. (A7) in the form 



2. Taylor Expansion for 



±^-fiR*^r{R*^r)), 



(A9) 



where /(i?*, r(_R*, t)) is the function at the right hand side of Eq. (A7). When one differentiates both sides of Eq. 
(A9) with respect to R* , one gets the result 



It is easy to show that 
df 



dR* 



= 2R* 



+ 



2M AAPiR*^ + 1) 
R*r^ 



1/2 



OR* dr dR* 



R*r 
2M 



4M2(i?*2 + 1) 



2M 4Af2(i?*2 ^ 1) 



-3i?*(i?*2 + 1)1/2 arccos 



2M 4M2(i?*2 ^ 

-1-1/2 

1/2- 



-1/2 



2M(i?*2 + 1) 



and 



a/ _ i?*2 + 1 

dr 2r 

i?*2 + 1 



2M 4Af2(i?*2 + 1) 

^2 



1/2 



AM 



2M 4M2(i?*2 + 1) 



-1/2 



Thus, according to Eq. (AlO), we get 

df/dR* 



dr 

OR* 



df/dr 



+12M^R* 



= 6MR* - 



R*r 

i?*2 + 1 



(i?*2 + 1) 1 



-I 1/2 



2Mr 4Af2 



arccos 



1/2- 



2A/(i?*2 + 1) 



Instead of obtaining the Taylor expansion directly from Eq. (A13), one may define 



(AlO) 



(All) 



(A12) 



(A13) 



w := R*\ 



(A14) 



(A15) 
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with simpler expressions. If one needs at most third -orde r terms of the Taylor expansion for ^fr, one has to calculate 
only the three first terms and the last term in Eq. ( A16 ). 
Putting w ~ and r — 2M one finds that 



AM. 



w=0 



(A17) 



When obtaining the next term for Eq. (A16) one has to use the chain rule: 

dg{w,T) dg{w,r) ^ dg{w,r) dr 

dw dw dr dw 

It should be noted that on the left hand side of Eq. ( A18| ) g is considered as a function of w and r, whereas on the 
right hand side g is considered as a function of w and r. It is straightforward to show that 



(A18) 



dg{w,r) 
dw 



+ 



3M 



3M f w + 1 1 



-1/2 



and 



dg{w,r) 
dr 



1 



\ 2Mr U'P 

3M 
r 



arccos 



1/2- 



2M{w + 1) 



(A19) 



w + 1 
3M{w + l) fw + 1 



1 



-1/2 



2Mr 4AP 



arccos 



1/2- 



2M{'w + l) 



(A20) 



When solving ^ one has to use a similar method that was used in solving Let us write Eq. (A7) in the form 

r 



(w,t)). 



When one differentiates both sides of this equation with respect to w, one gets 

^ dh dh dr 
dw dr dw 



(A21) 



(A22) 



It can be showed that 



dh 
dw 



2M 4M2(w + 1) 

^2 



1/2 



r 
4M 



8A'P{w + l) 



2M 4M2(it;-hl) 

-1/2 



-1/2 



2M AAP{w + l) 



+ + 1)^/^ arccos 



1/2-1 



2M{w + 1) 



(A23) 



and 



dl 
dr 



w + 1 
2r 



w ■ 



2M 4M2(w + 1) 

^2 



1/2 



4M 



2M 4M2(ui + l) 



-1/2 



(A24) 



which, together with Eq. (A22), lead to 

dr dh/dw 



dw 



dh/dr 



= 3M - 



+6]\/P 



1 



w + 1 



2Mr 4M2 



2(w + l) 

1/2 



1/2- 



2M(u; + 1) 



(A25) 
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Now, according to equations ([A18|), (|A19D, (|A20|) and (|A25D we get 



dw 



3r 



9M 9M2 



2(w + l)2 2(w + l) r 
18M3(w+ 1) 



■ arccos 



1/2- 



w + 1 I 



w + I \ 2Af r 4Af2 
36M3 / w + 1 I 



\ 2M(w + 1) 

1/2 

I arccos 



2M{w + 1) 



1/2 



r \ 2Mr 4A/2 



arccos 



1/2. 



1/2. 



2M{w + 1) 



Again, putting w — and r = 2Af one finds that 



dg 
dw 



= 3M. 



(A26) 



(A27) 



w=0 



To determine the next term in Eq. (A16) we must solve The chain rule, together with equations (A5) and 
(|A2C|), gives 



dg{w, t) dgiw, r) Or dg ( 2M 1 
dr dr dr dr \ r w + 1 



1/2 



2M 



6M2 



(u) + 1)3/2 + 1)1/2 

6M2(w + 1)1/2 



w + 1 1 



1/2 



arccos 



2Mr 4M2 
\ i/2n 



2M{w + 1) 



Thus, we get 



dr 



Let us now define 51 := A gain, using the chain rule and Eq. (|A5|) we get 



d'^g{w, t) dgi dgi dr dgi ( 2M 1 



1/2 



dr^ 



dr dr dr dr 



w + l 



One can show that 



dgi 
dr 



ZM{w + 1)1/2 



2r2(w + 1)1/2 



6M^ 



w + 1 1 



-1/2 



2Mr 4M2 



w + 1 1 



1/2 



r2(u; + 1)1/2 I 2Afr 471,^2 
12M2(u; + 1)1/2 



1/2- 



2M(u;+ 1) 



which, together with Eq. (ABO), leads to 

d^g{w,T) M 



dr^ 



+ 



6M2 



12M3 I w + 1 1 



(w + 1) r3 r'^{w + l)\2Mr 4M2 



24M3 (w + 1 1 



1/2 



2Mr 4M2 



1/2- 



2M{w + 1) 



(A28) 



(A29) 



(A30) 



(A31) 



(A32) 
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Putting again w = and r = 2M gives 



9t2 



ui=0 



1 



(A33) 



Finally, using equations ( A17D, ( A27D, (A29) and (A33) we may write down the Taylor expansion for Further- 



(A34) 



more, by using definitions ( Al4| ) and ( |A15 ) we may obtain the Taylor expansion for We get: 



dr 

OR* 



AMR* + 3Afi?* 



2M 



which is Eq. (p5|) 
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